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We have performed systematic large-scale all-electron correlated calculations on boron clusters 
B„(n=2-5), to study their linear optical absorption spectra. Several possible isomers of each clus- 
ter were considered, and their geometries were optimized at the coupled-cluster singles doubles 
(CCSD) level of theory. Using the optimized ground-state geometries, excited states of differ- 
ent clusters were computed using the multi-reference singles-doubles configuration-interaction 
(MRSDCI) approach, which includes electron correlation effects at a sophisticated level. These 
CI wave functions were used to compute the transition dipole matrix elements connecting the 
ground and various excited states of different clusters, eventually leading to their linear absorp- 
tion spectra. The convergence of our results with respect to the basis sets, and the size of the 
CI expansion were carefully examined. The contribution of configurations to many body wave- 
function of various excited states suggests that the excitations involved are collective, plasmonic 
type. 
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1. INTRODUCTION 

The area of cluster science has witnessed an enor- 
mous progress in terms of both the experimental 
and theoretical investigations of clusters of various 
atoms over last few decades. ^'^'^''^ Ranging from 
small clusters having a few atoms, to nanotubes, 
nanosheets, fullerenes etc., clusters have proven 
their usefulness in the fast emerging field of nan- 
otechnology.^ Because of the finiteness of the size, 
properties of clusters can be greatly tuned, and, 
therefore they are more amenable to nanoengineer- 



ing than their bulk counterpart.^ The evolution of 
clusters towards the bulk with the increasing num- 
ber of atoms, and the underlying mechanism, is a 
research topic of great contemporary interest.^' ^ 

At the present time, boron clusters are attract- 
ing great attention because of their novel proper- 
ties, and potential applications in nanotechnology 
and hydrogen storage related capabilities.^' ^"^'^^'^^ 
Boron atom, having s'^p^ valence electronic configu- 
ration, has short covalent radius and tends to form 
strong directional bonds producing clusters of cova- 
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lent nature. Because of this strong covalent bond- 
ing, it has hardness close to that of diamond. The 
ability of boron to form structures of any size due to 
catenation is only comparable to its neighbor car- 
bon.^ Planar boron clusters exhibit aromaticity^^ 
due to the presence of itinerant vr electrons, and 
some of them are analogous to aromatic hydrocar- 
bons.^'^ Boron fullerenes, boron sheets and single- 
sheet boron nitride — a graphene analogue — are the 
other examples of boron-based clusters. 

As far as the studies of boron-based clusters 
are concerned, small ionic boron clusters B^ (n ^ 
20) were experimentally studied by Hanley, Whit- 
ten and Anderson. Wang and coworkers have re- 
ported joint theoretical and experimental studies of 
the electronic structure of bare boron wheels, rings, 
tubes and large quasi-planar clusters.^^'^^'^^'^^ Us- 
ing the photoelectron spectroscopy, they predicted 
that tubular B20 can act as the smallest boron 
single walled nanotube. Transition metal-centered 
boron ionic ring clusters were studied by Constantin 
et. al.,^^ in a photo-electron spectroscopy exper- 
iment, supported by first-principles calculations. 
The abundance spectrum of boron clusters gener- 
ated by laser ablation of hexagonal boron nitride 
was studied by time of flight measurements per- 
formed by La Placa, Roland and Wynne. ^'^ They 
also postulated the existence of B36N24 cluster hav- 
ing a structure similar to that of Cgo fullerene. Lau- 
ret et al?^ probed the optical transitions in single 
walled boron nitride nanotubes by means of optical 
absorption spectroscopy. 

Larger pure boron clusters have also been in- 
vestigated extensively. Cage-like structure of Bgo — 
similar to Ceo fullerene — has been proposed theoret- 
ically.^^ A density functional theory (DFT) study of 
pure boron sheets and nanotubes was carried out by 
Cabria, Lopez and Alonso to explore their potential 
hydrogen storage materials. Chacko, Kanhere and 
Boustani investigated different equilibrium geome- 
tries of B24 cluster using Born-Oppenheimer molec- 
ular dynamics within the framework of DFT.^^ Ab- 
durahman et al.'^^ studied the ladder-like planar 
boron chains B„ (n=4-14), and computed their 
static dipole polarizabilities using the ab initio CI 
method. Johansson discussed strong toroidal ring 
currents in B20 and other toroidal boron clusters.^'' 
Double aromaticity was proposed in toroidal boron 
clusters B2n (n = 6,14) by Bean and Fowler. 

Regarding the smaller sized boron clusters, an 
early theoretical study of boron dimer was carried 



out by Langhoff and Bauschlicher,' who performed 
an extensive calculations using the complete-active- 
space self-consistent-field (CASSCF) multireference 
configuration interaction (MRCI) with a large ba- 
sis set. A similar study was carried out by Bruna 
and Wright^^ for the excited states of B2, and by 
Howard and Ray^® using the many-body perturba- 
tion theory. A systematic geometry and electronic 
structure calculations of bare boron clusters was 
reported by Boustani. He performed all-electron 
calculations at the SDCI level, but the contracted 
Gaussian basis sets used were small. Niu, Rao and 
Jena,^*^ using DFT and quantum chemical meth- 
ods, presented an account of electronic structures of 
neutral and charged boron clusters. In their study 
on small clusters, Moller-Plesset perturbation the- 
ory of fourth order (MP4) was used to account for 
the electron correlation effects. More recently, Ati§, 
Ozdogan, and Giiveng investigated structure and 
energetics of boron clusters using the DFT.'^^ Aro- 
maticity in planar boron clusters was addressed by 
Aihara, Kanno and Ishida.'^^ 

In spite of many theoretical studies of boron 
clusters of various shapes and sizes, very little ex- 
perimental information about their ground and ex- 
cited states is available. Conventional mass spec- 
trometry can distinguish between different clusters 
only according to their mass, but not according to 
their geometry. One has to rely on other theoreti- 
cal or experimental data to be able to differentiate 
one isomer from another. For example, using first 
principles calculations of vibrionic fine structure in 
C^g, and comparing it with experimentally avail- 
able data, Saito and Y. Miyamoto'^'^ identified the 
cage and bowl structures. Optical absorption spec- 
troscopy, coupled with extensive theoretical calcula- 
tions of the optical absorption spectra, can be used 
to distinguish between distinct isomers of clusters 
produced experimentally, because normally optical 
absorption spectra are sensitive to the geometries of 
the clusters. The optical absorption of alkali metal 
clusters has been extensively studied both experi- 
mentally and theoretically.^^' However, a 
very few such studies exist for the case of boron clus- 
ters. Marques and Botti'^^ calculated optical absorp- 
tion on different B20 isomers using time-dependent 
(TD) DFT. Boron fullerenes such as B38, B44, Bgo 
and B92 were also studied by Botti and coworkers^*^ 
using the same technique. However, to the best of 
out knowledge, there are no experimental and the- 
oretical study of optical absorption on other bare 
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boron clusters, particularly the smaller ones. It is 
with the aim of filling this void that we undertake a 
systematic study of the optical absorption in small 
boron clusters B„(n=2-5), employing the MRSDCI 
method, and high-quality Gaussian basis functions. 
We perform careful geometry optimization for each 
possible isomer, and compute the optical absorption 
spectra of various structures. We also analyze the 
many-body wave functions of various excited states 
contributing to the peaks in the computed spectra, 
and conclude that most of the excitations are collec- 
tive in nature, signalling the presence of plasmons. 

The remainder of this paper is organized as 
follows. Next section describes the theoretical and 
computational details of the work, followed by sec- 
tion 3 in which our results are presented and dis- 
cussed. In section 4 we present our conclusions and 
discuss possibilities for future work. Detailed infor- 
mation about various excited states contributing to 
optical absorption is presented in the Appendix. 

2. THEORETICAL AND 

COMPUTATIONAL DETAILS 

The geometry optimization of various isomers 
was done using the size-consistent coupled-cluster 
singles doubles (CCSD) based analytical gradi- 
ent approach, as implemented in the package 
GAMESS-us."^^ For the purpose, we used the 6- 
311G(d,p) basis set included in the program li- 
brary,^^ which is known to be well-suited for this 
task. The process of optimization was initiated by 
using the geometries reported by Ati§ et al.,"^^ based 
upon first principles DFT based calculations. For 
some simple geometries such as B2, B3 (D3/1 sym- 
metry), the optimization was carried out manually, 
by performing the MRSDCI calculations at difi^erent 
geometries, and locating the energy minima. Fig- 
ure 1 shows the final optimized geometries of the 
isomers studied in this paper. 

Once the ground state geometries of various 
isomers were determined, the correlated calcula- 
tions were performed to calculate their low-lying 
excited states using the multi-reference singles- 
doubles configuration-interaction (MRSDCI) ap- 
proach as implemented in the computer program 
MELD.''^ MRSDCI approach is a well-estabhshed 
quantum-chemical approach in which one consid- 
ers singly- and doubly-excited configurations from 
a number of reference configurations, leading to a 
good treatment of electron correlations, both for the 
ground, and the excited states, in the same calcu- 



lation. Using the ground- and excited-state wave 
functions obtained from the MRSDCI calculations, 
electric dipole matrix elements are computed and 
subsequently utilized to compute the linear absorp- 
tion spectrum assuming a Lorentzian line shape. By 
analyzing the wave functions of the excited states 
contributing to the peaks of the computed spectrum 
obtained from a given calculation, bigger MRSDCI 
calculations were performed by including a larger 
number of reference states. The choice of the ref- 
erence states to be included in a given calculation 
was based upon the magnitude of the corresponding 
coefficients in the CI wave function of the excited 
state (or states) contributing to a peak in the spec- 
trum. This procedure was repeated until the com- 
puted spectrum converged within an acceptable tol- 
erance, and all the configurations contributing sig- 
nificantly to various excited states were included in 
the list of the reference states. In the past, we have 
used such an iterative MRSDCI approach on a num- 
ber of conjugated polymers to perform large-scale 
correlated calculations of their linear and nonlinear 
optical properties.^'^'^^'^^''^^ 

The number of molecular orbitals, and thus the 
size of the CI expansion, increases rapidly with the 
increasing number of atoms in the clusters. Such 
a proliferation in the size of calculations can es- 
sentially make high-quality MRSDCI calculations 
impossible even for clusters of the sizes discussed 
in this work. Therefore, wherever possible, we have 
used the point-group symmetries corresponding to 
D2h, and its subgroups, at all levels of calculations 
to reduce the size of the CI expansions. During the 
MRSDCI calculatiosn, the frozen-core approxima- 
tion was employed, i.e., while constructing the CI 
expansion, no virtual excitations from the Is core 
electrons of the boron atoms of the cluster were con- 
sidered. Similarly, excitations into very high energy 
virtual orbitals were not considered with the pur- 
pose of keeping the calculations manageable. The 
impact of both the frozen-core approximation, and 
the deletion of high-energy virtual orbitals, along 
with the infiuence of the choice of the basis sets on 
our calculations will be examined carefully in the 
next section. 

3. RESULTS AND DISCUSSION 

In this section, first we discuss the convergence of 
our calculations with respect to various approx- 
imations and truncation schemes. Thereafter, we 
present and discuss the results of our calculations 
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Figure 1. (Color online) Geometry optimized structures of boron clusters with point group symmetry and the electronic 
ground state at the configuration interaction level. 

for varfous clusters. 

3.1. Convergence of Calculations 

Here, we carefully examine the convergence of the 
calculated absorption spectra with respect to the 
size and quality of the basis set, along with various 
truncation schemes in the CI calculations. 

3.1.1. Choice of the basis set 

In general, the results of electronic structure cal- 
culations depend upon the quality and the size of 
the basis set employed. While several contracted 
Gaussian basis functions have been devised which 
can deliver high-quality results on various quanti- 
ties such as the total energy, correlation energy, 
and the static polarizabilities of molecules, to the 
best of our knowledge the basis set dependence of 
linear optical absorption has not been explored. 
Since boron shows strong covalent bondings, the 
basis set used for calculations should have diffuse 
Gaussian contractions. Therefore, to explore the ba- 
sis set dependence of computed spectra, we used 
several basis sets^^'^*^ to compute the optical ab- 
sorption spectrum of the smallest cluster, i.e., B2. 
For the purpose, we used correlation-consistent ba- 
sis sets named AUG-CC-PVTZ, DAUG-CC-PVDZ, 
AUG-CC-PVDZ, CC-PVDZ, and DZP, which con- 
sist of polarization functions along with diffuse ex- 



ponents, and were designed specifically for post 
Hartree-Fock correlation calculations.^^' From the 
calculated spectra presented in Fig. 2 the follow- 
ing trends emerge: the spectra computed by vari- 
ous augmented basis sets (AUG-CC-PVTZ, DAUG- 
CC-PVDZ, AUG-CC-PVDZ) are in good agree- 
ment with each other in the energy range up to 8 
eV, while those obtained using the nonaugmented 
sets (CC-PVDZ and DZP) disagree with them sub- 
stantially, particularly in the higher energy range. 
Given the fact that augmented basis sets are consid- 
ered superior for molecular calculations, we decided 
to perform calculations on the all the clusters using 
the AUG-CC-PVDZ basis set. This is the smallest 
of the augmented basis sets considered by us, and, 
therefore, does not cause excessive computational 
burden when used for larger clusters. 
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Figure 2. (Color online) Optical absorption in B2 calculated 
using various Gaussian contracted basis sets. Increasing more 
and more diffuse d type Gaussians shows negligible effect on 
optical spectra. 

3.1.2. Orbital Truncation Schemes 

If the total number of orbitals used in a CI expan- 
sion is N, the number of configurations in the cal- 
culation proliferates as ~ N^, which can become 
intractable for large values of A^. Therefore, it is 
very important to reduce the number of orbitals 
used in the CI calculations. The occupied orbitals 
are reduced by employing the so-called "frozen-core 
approximation" described earlier, while the unoccu- 
pied (virtual) set is reduced by removing very high- 
energy orbitals. 



CORE ACTIVE 
CORE FROZEN 




Figure 3. (Color online) The effect of freezing the core or- 
bitals (Is) of boron atoms on optical absorption spectrum of 
B2. It renders almost no effect on optical absorption spec- 
trum. 

The influence of freezing the Is core orbitals 
on the optical absorption spectrum of B2 cluster is 
displayed in Fig. 3, from which it is obvious that it 
makes virtually no difference to the results whether 
or not the core orbitals are frozen. The effect of re- 
moving the high-energy virtual orbitals on the ab- 



sorption spectrum of B2 is examined in Fig. 4. From 
the figure it is obvious that if all the orbitals above 
the energy of 1 Hartree are removed, the absorption 
spectrum stays unaffected. Therefore, in rest of the 
calculations, wherever needed, orbitals above this 
energy cutoff were removed from the list of active 
orbitals. Theoretically speaking this cutoff is sound, 
because we are looking for absorption features in the 
energy range much smaller than 1 Hartree. 




Figure 4. (Color online) The effect of the number of active 
orbitals (Nact) on the optical absorption spectrum of B2. 
Until Nact=42, the optical spectrum does not exhibit any 
significant change. It corresponds to 1.0 Hartree (~ 27.2 eV) 
virtual orbital energy. 



3.1.3. Size of the CI expansion 

As mentioned earlier that the electron correlation 
effects in both the ground and the excited states 
were accounted in our calculations by including the 
relevant configurations in the reference list of the 
MRSDCI expansion. The greater numerical accu- 
racy demands the inclusion of a large number of 
configurations in the reference list, but that leads to 
a rapid growth in the size of the CI expansion, mak- 
ing the calculations numerically prohibitive. How- 
ever, here we are interested in computing the en- 
ergy differences rather than the absolute energies 
of various states, for which good accuracy can be 
achieved even with moderately large CI expansions. 
In Table 1 we present the average number of refer- 
ence states (Nj-ej) included in the MRSDCI expan- 
sion and average number of configurations (N total) 
for different isomers. For a given isomer, the aver- 
age has been calculated across different irreducible 
representations which were needed in these symme- 
try adapted calculations in order to compute the 
ground and various excited states. The extensive- 
ness of our calculations can be seen from the number 
Ntotah which is « 77000 for the simplest cluster, and 
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around four million for each symmetry subspace of 
B5. This makes us believe that our results are fairly 
accurate. 

Before we discuss the absorption spectrum for 
each isomer, we present the ground state energies 
along with the relative energies of each isomer are 
given in Table 2. The MRSDCI energy convergence 
threshold was 10"^ for all the isomers, with 10 
as convegrence threshold for configuration coeffi- 
cients. From the results it is obvious that as far 
as the energetics are concerned, for the B3 the tri- 
angular structure is most stable, while for B4 and 
B5 the rhombus and pentagonal structures, respec- 
tively, are favorable. 

Table 1. The average number of reference con- 
figurations (Nj.gy), and average number of total 
configurations total) involved in MRSDCI cal- 



culations 


of various isomers of 


boron 


clusters. 


Cluster 


Isomer 




^total 


B2 


Linear 


24 


77245 


B3 


Triangular 


36 


596798 




Linear 


41 


671334 


B4 


Rhombus 


37 


1127918 




Square 


40 


1070380 




Linear 


34 


1232803 




Distorted Tetrahedron 


28 


1253346 


B5 


Pentagon 


22 


3936612 




Distorted Tri. bipyramid" 


7 


3927508 



Note: ° Cs symmetry of isomer converted to Ci in calcu- 
lations. 



Table 2. Ground state (GS) energies (in Hartree) at MRSDCI 
level and the relative energies (in eV) of different isomers of clus- 
ters. 



Cluster 


Isomer 


GS energy 


Relative 






(Ha) 


energy (eV) 


B2 


Linear 


-49.27844 


0.00 


B3 


Triangular 


-73.98998 


0.00 




Linear 


-73.92906 


1.66 


B4 


Rhombus 


-98.74004 


0.00 




Square 


-98.73785 


0.06 




Linear 


-98.66575 


2.02 




Distorted Tetrahedron 


-98.63213 


2.94 


B5 


Pentagon 


-123.42652 


0.00 




Distorted Tri. Bipyramid 


-123.31485 


3.04 



3.2. Calculated Photoabsorption 
Spectra of Various Clusters 

Next we present and discuss the results of our pho- 
toabsorption calculations for each isomer. 



3.2.1. B2 

The simplest and most widely studied cluster of 
boron is B2 with Doofc point group symmetry. Us- 
ing the SDCI method, we obtained its optimized 
bond length to be 1.59 A (c/. Fig. 1(a)), which is 
in excellent agreement with the experimental value 
1.589 A..'''^ Using a DFT based methodology, Ati§ 
et al.,^^ obtained a bond length of 1.571 A, while 
Howard and Ray calculated it to be 1.61 A, using 
the fourth-order perturbation theory (MP4).^^ 

Because the ground state of B2 is a spin triplet, 
its many-particle wave function predominantly con- 
sists of a configuration with two degenerate singly- 
occupied molecular orbitals (SOMO) referred to as 
Hi and H2 in rest of the discussion. The excited 
state wave functions will naturally consist of con- 
figurations involving electronic excitations from the 
occupied MOs to the unoccupied MOs starting from 
lowest unoccupied molecular orbital (LUMOs, L for 
short). Our calculated photoabsorption spectrum 
shown in Fig. 5 is characterized by weaker absorp- 
tions at low energies, and a very intense one at high 
energy. The many-particle wave functions of excited 
states contributing to various peaks are presented in 
Table 3. A feeble peak appears near 0.85 eV, domi- 
nated by H2 — )• L and H2 L + A excitations com- 
pared to the HF reference configuration. It is fol- 
lowed by a couple of smaller peaks at 4.20 and 4.91 
eV. The most intense peak is found at 7.05 eV, to 
which two closely spaced states contribute. Transi- 
tion to the state near 6.97 eV is polarized transverse 
to the bond length, while the one close to 7.05 eV 
carries the bulk of oscillator strength, and is reached 
by longitudinally polarized photons. All these states 
exhibit strong mixing of singly-excited configura- 
tions. Near 8 eV, a smaller peak appears which has 
strong contributions from doubly-excited configu- 
rations H — I ^ L; Hi ^ L + 2 and — 1 — )• L; 
H2 — )• L+2. The wave functions of the excited states 
contributing to all the peaks exhibit strong configu- 
ration mixing, instead of being dominated by single 
configurations, pointing to the plasmonic nature of 
the optical excitations.^^ 
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Figure 5. The linear optical absorption spectrum of B2, cal- 
culated using the MRSDCI approach. The peaks correspond- 
ing to the light polarized along the molecular axis are labeled 
with the subscript ||, while those polarized perpendicular to 
it are denoted by the subscript _L. For plotting the spectrum, 
a uniform linewidth of 0.1 eV was used. 

3.2.2. B3 

Boron trimer has two possible isomers, triangular 
and the linear one shown in Figs. 1(b) and 1(c). We 
found equilateral triangle with D^h symmetry to be 
the most stable isomer. The optimized bond length 
for triangular isomer is 1.55 A, with the ground 
state (^^1) energy 1.66 eV lower than that of its 
linear counterpart. We also explored the possibility 
of isosceles triangular structure favorable one, 
because B3 is an open-shell system, making it a pos- 
sible candidate for Jann- Teller distortion. However, 
the CCSD optimised geometry corresponding to the 
isosceles structure is so slightly different compared 
to the equilateral one, that it is unlikely to affect the 
optical absorption spectrum in a significant manner. 
Our calculated bond length is in good agreement 
with experimental value 1.57 A,^^ as well as with 
other reported theoretical values of 1.553 A,^^ 1.56 
A28 and 1.548 A.^i 

The linear B3 isomer with the T)ooh symme- 
try, and the as ground state, was found to 
have equal bond lengths. Our SDCI optimized bond 
length of 1.51 A agrees well with the value 1.518 A 
reported by Ati§ et al.^^ 

The photoabsorption spectra of two isomers of 
B3 are presented in Figs. 6 and 7. The correspond- 
ing many-particle wave functions of excited states 
contributing to various peaks are presented in Ta- 
ble 5 and 4. It is obvious that in the linear struc- 
ture, absorption begins at a lower energy as com- 
pared to the triangular one, although the intensity 



of its low-energy peaks is very small. In the trian- 
gular isomer on the other hand, most of the inten- 
sity is concentrated at rather high energies, except 
for a weaker peak close to 3 eV. The optical spec- 
tra of linear isomer begins with very weak peaks 
at 0.7 eV (longitudinal polarization) and 2.7 eV 
(transverse polarization), with their many-particle 
wave functions dominated by singly-excited config- 
urations. The relatively intense peak at 4.3 eV cor- 
responding to a longitudinally polarized transition, 
is dominated by doubly-excited configuraitons. It is 
followed by a small peak mainly due to single excita- 
tion H — 2 —7- L, near 5.9 eV. The most intense peak 
of the spectrum occurs at 7.4 eV, followed by an- 
other strong peak close to 7.7 eV. Both the features 
correspond to longitudinally polarized transitions, 
with the many particle wave functions of the con- 
cerned states being strong mixtures of single and 
double excitations with respect to the HF reference 
state. We note that in the absorption spectrum of 
the linear cluster, quite expectedly, the bulk of the 
oscillator strength is carried by longitudinally po- 
larized transitions. 

Because the triangular cluster is a planar clus- 
ter, its orbitals can be classified as in-plane a or- 
bitals, and the out-of-plane vr orbitals. Both the 
HOMO (a singly occupied orbital, in this case) and 
the LUMO for this isomer are cr-type orbitals. For 
this system, two types of optical absorptions are 
possible: (a) those polarized in the plane of the clus- 
ter, and (b) the ones polarized perpendicular to that 
plane. Our calculations reveal that the transitions 
correponding to perpendicular polarization {z direc- 
tion), except for a couple of peaks, have negligible 
intensities. From Fig. 7 it is obvious that the opti- 
cal absorption in the triangular isomer starts with 
a very weak z— polarized feature near 0.8 eV (peak 
I), corresponding to a state with the wave function 
dominated by single excitations (vr — )• a*). This is 
followed by a series of peaks ranging from II to VI 
which correpsond to the photons polarized in the 
plane of the cluster. All these peaks are dominated 
by states consisting primarily of singly-excited con- 
figurations of the (T — )• cr* type. The most intense 
peak VI is followed by a shoulder-like feature (VII) 
corresponding to a z-polarized absorption. 

If we compare the absorption spectra of the lin- 
ear and the triangular B3, the peak at 4.34 eV in the 
spectrum of the linear cluster is the distinguishing 
feature, and can be used to differentiate between 
the two isomers. 
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Figure 6. The linear optical absorption spectrum of linear 
B3, calculated using the MRSDCI approach. The peaks cor- 
responding to the light polarized along the molecular axis 
are labeled with subscript ||, while those polarized perpen- 
dicular to it are denoted by the subscript _L. For plotting the 
spectrum, a uniform linewidth of 0.1 eV was used. 




Figure 7. (Color online) The linear optical absorption spec- 
trum of triangular B3 calculated using the MRSDCI ap- 
proach. Peaks corresponding to light polarized in the plane of 
the molecule are labeled with subscript ||, while those polar- 
ized perpendicular to the plane are denoted by the subscript 
_L. For plotting the spectrum, a uniform linewidth of 0.1 eV 
was used. 



3.2.3. B4 

For the B4 cluster, we investigated the rhombus, 
square, linear and tetrahedral structures. While the 
rhombus shaped isomer was found to have the low- 
est energy, but the square isomer is higher in en- 
ergy only by a small amount. As a matter of fact, 
at the HF level the energies of the two isomers were 
found to be almost degenerate. It was only after the 
electron correlation effects were included at the CI 
level that the rhombus stabilized by ~0.06 eV (c/. 
Table 2)with respect to the square. For the rhom- 



bus, the ground state was ^Ag, with the optimized 
bond length 1.529 A, and the short diagonal length 
1.909 A. These results are in good agreement with 
with the correponding lengths of 1.528 A and 1.885 
A reported by Boustani,,^^ and 1.523 A and 1.884 A 
computed by Ati§ et al.^^ Both HOMO and LUMO 
of rhombus isomer are a orbit als. 

For the square isomer, with D4/1 symmetry, 
the electronic ground state is expectedly ^Aig. As 
shown in Fig. 1(e), our optimized bound length is 
1.53 A, which agrees well with the values 1.527 A 
and 1.518 A as reported in Refs. 29 and 31. In this 
isomer, HOMO is a cr orbital while LUMO is a tt 
orbital. 

Linear B4, with the T)ooh symmetry, has the 
electronic ground state of ^S". However, energet- 
ically linear structure is 2.02 eV higher than the 
rhombus one (c/. Table 2) which rules out its ex- 
istence at the room temperatures. As per Fig. 1(f) 
, the central bond length was found to be 1.49 A, 
with the two outer bonds being 1.55 A in length. For 
the same bonds, Ati§ et al. reported these lengths 
to be 1.487 A and 1.568 A, respectively.^^ 

The distorted tetrahedral structure having Cs^ 
symmetry, made up of four isosceles triangular faces 
with lenghts 1.785 A, 1.785 A and 1.512 A. This iso- 
mer also lies much higher in energy as compared to 
the most stable rhombus structure. 

The absorption spectra of rhombus, square, lin- 
ear, and tetrahedral isomers are presented in Figs. 
8, 9, 10, and 11 respectively. From the figures it is 
obvious that the general features of the absorption 
spectra of rhombus and square isomers are similar, 
except that the rhombus spectrum, with the on- 
set of the absorption near 4 eV, is red-shifted by 
about 1 eV as compared to the square. The ab- 
sorption spectrum of the linear structure is slightly 
red-shifted as compared to the rhombus and square 
shaped isomers, with the majority of absorption oc- 
curing in the energy range 5-8 eV. This aspect of 
the photoabsorption in B4 is similar to the case of 
B3 for which also the linear structure exhibited a 
redshifted absorption compared to the triangular 
one. 

Since B4 rhombus isomer has D2h symmetry, we 
can represent the absorption due to light polarized 
in different directions in terms of irreducible rep- 
resentations of D2h- So absorption due to in-plane 
polarized light corresponds to Biu and i?2?i, while 
B'iu corresponds to light polarized in the direction 
perpendicular to the plane of the isomer. 
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The polarization resolved absorption spectrum 
of rhombus B4, as shown in Fig. 8, exhibits a rather 
blue-shifted nature as compared to the linear iso- 
mer. The many-particle wave functions of excited 
states contributing to various peaks are presented 
in Table 6. The onset of spectrum is seen at 4.15 
eV followed by a peak at around 6.12 eV. Both of 
them are due to y— polarized component, i.e. along 
the larger diagonal. The dominant contribution to 
these peaks come from o" — )• vr* for former, and 
TT — )• TT* for latter. The x-component does not con- 
tribute much in the whole spectrum, except for mi- 
nor peaks at 4.2 eV and 6.6 eV. It is characterised 
by mainly vr — )• vr* type transitions. It is followed 
by a relatively low intensity peak at 7.3 eV due 
to y-polarized component with leading contribution 
from 0" — )• vr* transitions. The most intense peak, at 
7.84 eV, having y-polarization component, is char- 
acterized by o" — )• vr* type of transitions. There are 
no direct H ^ L transitions for this isomer, be- 
cause they are dipole forbidden. The absorption due 
to light polarized in the direction perpendicular to 
the plane of isomer is negligible. 
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Figure 8. (Color online) The linear optical absorption spec- 
trum of B4 rhombus geometry using the MRSDCI approach. 
Isomer is aligned in x — y plane with short diagonal along 
2;-axis. Peaks corresponding to light polarized along x and 
y-axis are labeled with subscript x and y. For plotting the 
spectrum, a uniform linewidth of 0.1 eV was used. 

The square B4 isomer, because of its symme- 
try, gets equal contribution to absorption spectrum 
from both x— and y— component. It corresponds to 
in-plane polarization due to Biu and i?2u irreducible 
representation, while B^u corresponds to light po- 
larized in the direction perpendicular to the plane 
of the isomer. However, in this isomer also, the 
contribution due to latter is quite negligible. The 
many-particle wave functions of excited states con- 



tributing to various peaks are presented in Table 
7. It shows just one major peak at 4.88 eV below 
7 eV, characterized by o" — t- 7r*;o" — t- vr* double ex- 
citation. Two smaller peaks appear in this range 
at 5.5 eV and 6.4 eV, with leading contributions 
from o" — )• 7r*;cj — 7- vr* and cr — >• vr*;vr — t- vr* exci- 
tations respectively. Beyond 7 eV, there are many 
closely spaced peaks including the most intense one 
at 7.89 eV. It is characterized by double excita- 
tion a —7- vr*;vr — t- vr*. In this isomer also, a direct 
H ^ L transitioin is forbidden. Though, there is 
very little difference in total energies of rhombus 
and square isomers of B4, their optical absorption 
spectra are completely different. They can be eas- 
ily identified from each other by looking at number 
of peaks below 7 eV energy. Rhombus exhibits two 
major peaks, while square has just one. 




E (eV) 

Figure 9. (Color online) The linear optical absorption spec- 
trum of B4 square geometry using MRSDCI approach. Iso- 
mer is aligned in x — y plane. Spectrum represents the equal 
contribution from light polarized in x and y direction. Peaks 
corresponding to light polarized in the plane of the molecule 
are labeled with subscript ||. For plotting the spectrum, a 
uniform linewidth of 0.1 eV was used. 



Linear B4 isomer exhibits absorption with few, 
but sharp peaks. The many-particle wave functions 
of excited states contributing to various peaks are 
presented in Table 8. The onset of optical absorp- 
tion occurs near 4.5 eV, due to absorption of long- 
axis polarized light, followed by two major peaks 
at 5.95 eV, and 7.36 eV. The first of these two in- 
tense peaks, peak II is dominated by singly-excited 
configurations, while the second one (peak III) is a 
strong mixture of both singly- and doubly-excited 
configurations with respect to the HF reference con- 
figuration. 
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Figure 10. The optical absorption spectrum of linear B4, 
calculated using the MRSDCI approach. The peaks corre- 
sponding to the light polarized along the molecular axis are 
labeled with subscript ||, while those polarized perpendicu- 
lar to it are denoted by the subscript _L. For plotting the 
spectrum, a uniform linewidth of 0.1 eV was used. 

The 3D structure, a distorted tetrahedral iso- 
mer, exhibits an absorption spectrum very different 
from other isomers, as displayed in Fig. 11. The 
many-particle wave functions of excited states con- 
tributing to various peaks are presented in Table 
9. It is the only B4 isomer to exhibit peaks be- 
low 4 eV. The absorption spectrum is spread over 
a much larger energy range, and is almost contin- 
uous. The oscillator strengths associated with vari- 
ous peaks are much smaller than in other isomers, 
and most of the peaks appear pairwise. The onset 
of absorption spectrum is seen at around 1.1 eV, 
characterised mainly by an excited dominated by 
single-excitation H ^ L (cf. Table 9). In this iso- 
mer, in contrast to other B4 isomers, direct H ^ L 
transitions are allowed. Higher energy peaks in this 
isomer are dominated by doubly-excited configura- 
tions, and, are, therefore, sensitive to the electron- 
correlation effects. 
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Figure 11. (Color online) The linear optical absorption 
spectrum of B4 distorted tetrahedral geometry using the 
MRSDCI approach. Peaks corresponding to light polarized 
along X, y and 2:-axis are labeled with subscript x, y and z 
respectively. For plotting the spectrum, a uniform linewidth 
of 0.1 eV was used. 



3.2.4. ^5 

We investigated two isomers of B5: a Jahn- Teller 
distorted pentagon with the C2V symmetry, and (b) 
a triangular bipyramid with the point group 
symmetry. The latter one is the second 3-D struc- 
ture of the boron clusters probed in this work. The 
lowest lying pentagon isomer, has electronic 
ground state, and is 3.04 eV lower in energy as 
compared to the bipyramid structure. For the pen- 
tagon, the symmetry of ground state at the SCF 
level was Ai, however, at the MRSDCI level the B2 
state became lower in energy, in agreement with the 
previous calculations of Boustani.^^ Our optimized 
geometry for the pentagon (Fig. 1(h)) corresponds 
to an average bond length of 1.56 A, as against 1.57 
A reported by Boustani,^^ and 1.644 A reported by 
Ati§ et al.^^ The singly occupied molecular orbital 
(denoted by H) and LUMO of pentagon isomers are 
of vr and a type, respectiely. The bond lengths for 
the bipyramid structure are shown in Fig. l(i), with 
an average bond length of 1.704 A. The triangular 
base was found to be isosceles with 1.97 A as equal 
sides, and 1.75 A as the other side. 

The absorption spectra of the two isomers are 
presented in Figs. 12 and 13. The many-particle 
wave functions of excited states contributing to var- 
ious peaks are presented in Table 10 and 12 respec- 
tively. From the figures it is obvious that the absorp- 
tion in the bipyramid starts at much lower energies 
as compared to the pentagonal isomer. Intense ab- 
sorption peaks in pentagon B5 are located at en- 
ergies higher than 5 eV, with three equally intense 
peaks at 5.58 eV, 6.30 eV and 7.16 eV, with the 
photons polarized along the plane of the molecule 
direction. It has an underlying low intensity absorp- 
tion contribution from photons polarized along z— 
direction, which is perpendicular to the molecular 
plane. The major contribution to the peak at 5.58 
eV comes from vr — )■ cr* and cr — )• cr*, single exci- 
tations. The latter configuration also contributes to 
the most intense peak at 6.30 eV. The peak at 7.16 
eV is mainly due to ct — )• a* type transitions. 
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Figure 12. (Color online) The linear photo-absorption spec- 
trum of pentagon B5, calculated using the MRSDCl ap- 
proach. The peaks corresponding to the light polarized in 
the plane of the molecule are labeled with subscript ||, while 
those polarized perpendicular to it are denoted by the sub- 
script _L. For plotting the spectrum, a uniform linewidth of 
0.1 eV was used. 



Since the B5 triagonal bipyramid isomer is not 
a symmetric one, the calculations were done using 
Ci symmetry, thereby increasing the difficulty in 
diagonalizing the Hamiltonian. Hence, in order to 
reduce the matrix size, we have used a smaller num- 
ber of reference configurations, and also relaxed the 
energy convergence threshold criterion a little. 

The optical absorption spectrum of B5 triangu- 
lar bipyramid isomer is exhibited by almost equally 
spaced peaks at relatively lower energies. The opti- 
cal absorption starts at 1.74 eV characterized by 
H — 1 L + A configuration. It is followed by 
two equal intensity peaks at 3.16 eV and 4.27 eV 
with contributions from single excitations i7 — 2 — )• 
H and H—\ — )• L-l-2, respectively. The most intense 
peak is found at 7.52 eV dominated by the doubly 
excited configuration H — 2 ^ H ; H — 1 ^ L + 2. 
There are two distinguishing features as far as the 
optical absorption in the two isomers is concerned: 
(a) presence of three intense peaks in the higher en- 
ergy region of the absorption spectrum of the pen- 
tagonal isomer, and (b) occurance of equally spaced 
absorption peaks at lower energies in the spectrum 
of bipyramidal isomer. 
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Figure 13. (Color online) The linear photo-absorption spec- 
trum of distorted triangular bipyramid B5, calculated using 
the MRSDCl approach. Peaks corresponding to light polar- 
ized along X, y and z-axis are labeled with subscript x, y and 
z respectively. For plotting the spectrum, a uniform linewidth 
of 0.1 eV was used. 

4. Conclusions and Outlook 

We presented systematic large-scale all-electron cor- 
related calculations of photoabsorption spectra of 
boron clusters B„, (n =2-5) with several possi- 
ble isomers of each cluster. The calculations were 
perfomred using the MRSDCl method which takes 
electron correlations into account at a sophisti- 
cated level, both for the ground and the excited 
states. For a cluster consisting of a given number 
of atoms, significant changes were observed in ab- 
sorption spectra for different isomers, indicating a 
strong structure-property relationship. Therefore, 
our computed spectra can be used in the future pho- 
toabsorption experiments to distinguish between 
different isomers of a cluster, something which is 
not possible with the conventional mass spectrome- 
try. We also analyzed the many-particle wave func- 
tions of various excited states and found them to be 
a mixture of a large number of configurations, indi- 
cating the nature of photoexcited states in these 
clusters to be plasmonic.^^ A noteworthy aspect 
of the groundstate photoabsorption of various clus- 
ters was the absence of high-intensity peaks in the 
low-energy region of the spectrum. The most in- 
tense peaks occurred at higher energies involving 
orbitals away from the Fermi level, consistent with 
the fact that the bulk boron is an indirect bandgap 
semiconductor, with no optical absorption at the 
gap. In other words, optical absorption features of 
bulk boron were already evident in smaller clusters. 
Given the fact that aluminum has the same valence 
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shell structure as boron, it will be interesting to 
perform a similar set of calculations on small alu- 
minum clusters, in order to compare and contrast 
their photophysics. Calculations along those direc- 
tions are presently underway in our group, and re- 
sults will be communicated in future publications. 
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Appendix: Excited State CI Wavefunc- 
tions, Energies and Transition Moments 

In the following tables we have given the exci- 
tation energies (with respect to the ground state), 
the many body wavefunction and the oscillator 
strengths of the excited states corresponding to the 
peaks in the photoabsorption spectra of various iso- 
mers listed in Fig. 1, and discussed in section 3. 
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Table 3. Excitation energies, E, and many-particle wave functions of ex- 
cited states corresponding to the peaks in the linear absorption spectrum 
of B2 (c/. Fig. 5), along with the oscillator strength of the transitions 
(/12 = §^(£'2 — K'^Mi|G')P)i where, |m) denotes the excited state in 

question, \G), the ground state, and dj is the i-th Cartesian component of the elec- 
tric dipole operator. Longitudinal and transverse polarization corresponds to the 
absorption due to light polarized along and perpendicular to the molecular axis 
respectively. In the wave function, the bracketed numbers are the CI coefficients of 
a given electronic configuration. Symbols Hi,H2 denote SOMOs discussed earlier, 
and H, and L, denote HOMO and LUMO orbitals respectively. HF denotes the 
Hartree-Fock configuration. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 












\HF) (0.8673) 

\H -1 ^ L;H -l^ i)(0.2706) 
\H - 1 ^ L; H - I ^ L + 4)(0.1283) 


I 


0.845 


0.1187 


transverse 


\H2 L) (0.8742) 
\H2 ^ L + 4) (0.2194) 




II 


4.207 


1.1449 


longitudinal 


\H~1^ L)(0.7599) 
\Hi ^ L 4- 7) (0.2873) 
\H2 ^ L + 7) (0.2873) 




III 


4.914 


0.4252 


transverse 


\H-1^ L;H - 1 ^ 
Ti"- 2 ^i + 2) (0.3186) 
\H- 1 ^ L + 7)(0.266) 


-2) (0.7553) 


IV 


6.97 
7.05 


4.3060 
18.0923 


transverse 
longitudinal 


\H2 ^ L-H 4) (0.5600) 
\H2 L-H 5) (0.5150) 
\H-1^ L) (0.4462) 
\H2 -> L 4- 3) (0.3346) 
\Hi ^ L-h 3) (0.3346) 
1 ^ L + 4)(0.2921) 




V 


7.973 


0.6907 


longitudinal 


\H-l^L;H2 ^ L-h 2) (0.5613) 
\H-1^ L;Hi -> L-^- 2) (0.5613) 
\H- 1 ^ L + 4)(0.2303) 



Note: °GS does not correspond to any peak, instead it corresponds to the ground 
state wavefunction of the isomer. 
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Table 4. Excitation energies, E, and many-particle wave functions of ex- 
cited states corresponding to the peaks in the linear absorption spectrum 
of triangular B3 (c/. Fig. 7), along with oscillator strength of transition 
(/12 = §^(-^2 — El) |(m|di|G')p), where, |m) denotes the excited state 
in question, \G), the ground state, and dj is the i-th Cartesian component of 
the electric dipole operator. The polarization || corresponds to the absorp- 
tion due to Hght polarized in the plane of isomer, while _L corresponds to the 
polarization perpendicular to that plane. In the wave function, the bracketed 
numbers are the CI coefficients of a given electronic configuration. Symbol 
L and H denote LUMO and SOMO orbitals discussed earlier. HF denotes 
the Hartree-Fock configuration. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 










\HF) (0.8229) 
|i/-3^L) (0.2162) 
\H L)(0.2072) 
|i/ — 3 — > f/)(0.1544) 


I 


0.797 


0.0912 


_L 


Iff - 1 ^ Jf) (0.7357) 
\H -1^ L) (0.4693) 
\H -1 ^ L + 8)(0.1142) 


II 


1.710 


0.0266 
0.0282 


" 


\H ^ L-Hl)(0.4816) 
if- 2 ^Jf) (0.4315) 
ff-2^L) (0.3757) 
|if-3^L4-l) (0.3709) 
\H L)(0.5041) 
if — 3 ^ if) (0.4703) 


III 


2.840 
2.872 


0.5173 
0.4741 




|ii- 2^ if) (0.5374) 
|ii ^ L + 1)(0.5199) 
ii- 2 ^L) (0.3306) 
|ii-3^L + l)(0.2216) 
|ii ^ L)(0.4906) 
|ii-3^L) (0.4499) 


IV 


5.710 
5.730 


0.3259 
0.2053 


II 


if — 2 — )> L)(0.4114) 
|ii- 2^ if) (0.34 74) 
1 ff -3 ^L) (0.4202) 
ff-2->if;if-3-> L)(0.2050) 


V 


5.988 
6.02 


0.5051 
0.9679 




lii-3->L-hl)(0.3968) 
1 ii- 2 ^L) (0.3057) 
|ii- 2^ if) (0.1578) 
1 ff -3 ^L) (0.2578) 
Iff ^ L)(0.1908) 


VI 


7.657 
7.697 


1.6393 
1.0733 




Iff ^ L + 2) (0.5863) 
Iff -3 ^L-h 2) (0.2497) 
Iff ^ L + 14)(0.2431) 
Iff ^ L + 4) (0.2093) 
Iff ^ L + 3)(0.2991) 
1 ff -3 ^L + 1) (0.2648) 
1 ff-2^ if) (0.1794) 
ff-2^L) (0.1521) 


VII 


7.893 


2.1315 


_L 


Iff ^ L + 5)(0.7611) 
1 ff -3 ^L + 4) (0.3477) 



Note: "GS does not correspond to any peak, instead it corresponds to the 
ground state wavefunction of the isomer. 
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Table 5. Excitation energies, E, and many-particle wave functions of 
excited states corresponding to the peaks in the linear absorption spec- 
trum of linear B3 (c/. Fig. 6), along with oscillator strength of transition 
(/12 = §^(-£'2 — Ei)Y2i\{^\di\G)\'^), where, |m) denotes the excited state 
in question, |G), the ground state, and dj is the i-th Cartesian component of the 
electric dipole operator. Longitudinal and transverse polarization corresponds to 
the absorption due to light polarized along and perpendicular to the molecular 
axis respectively. In the wave function, the bracketed numbers are the CI co- 
efficients of a given electronic configuration. Symbols H and L denote HOMO 
and LUMO orbitals respectively, and Hi denotes SOMOs discussed earlier. HF 
denotes the Hartree-Fock configuration. 



Peak 


E (eV) 




Polarization 


Wave Function 




GS" 








\HF) (0.6650) 












\H - 


1 - 


L-H -1^ 


L) (0.3286) 










\H - 


1 - 


L-H -1^ 


L) (0.3277) 










\H - 


> L 


Li 


- 1)(0.2158) 










\H - 


> L;H-1^ Li 


- 1)(0.2157) 


I 


0.723 


0.1520 


longitudinal 


\H - 


> L) (0.8222) 




II 


2.707 


0.1337 


transverse 


\H - 


1 ^ 


L) (0.5360) 




III 


4.338 


4.5187 


longitudinal 


\H - 


1 - 


L-H -1^ 


i?) (0.5826) 










\H - 


1 - 


L;H -1^ 


if) (0.5826) 


IV 


5.937 


0.2295 


longitudinal 


\H - 


2 ^ 


L) (0.2497) 




V 


7.359 


8.4221 


longitudinal 


\H - 


1 - 


L + 1) (0.3683) 










\H - 


1 - 


L + 1) (0.3683) 










\H - 


3 


L-H ^ i)(0.3301) 










\H - 


3 - 


Jf) (0.2364) 




VI 


7.731 


2.9881 


longitudinal 


\H - 


2 - 


L) (0.3858) 












\H - 


> L 


;ii'-l^L-|^4) (0.2851) 










\H - 


> L-H-l ^ L-f 4) (0.2851) 




7.786 


2.3331 


transverse 


\H - 


1 ' 


L;H -1^ 


L + 4) (0.6438) 










\H - 


1 - 


L;H -1^ 


L + 4) (0.3905) 



Note: °GS does not correspond to any peak, instead it corresponds to the ground 
state wavefunction of the isomer. 



16 Ravindra Shtnde and Alok Shukla 



Table 6. Excitation energies, E, and many-particle wave lunctions 
ol excited states corresponding to the peaks in the Hnear absorption 
spectrum ol rhombus B4 (c/. Fig. 8), along with oscillator strength ol 
transition (/12 = § ^ (-E2 — -Ei) ^I] J {?Ti|cii|G) |^), where, |m) denotes 
the excited state in question, \G), the ground state, and is the 
i-th Cartesian component ol the electric dipole operator. The polar- 
ization in-plane corresponds to the absorption due to light polarized 
in the plane ol isomer. In the wave lunction, the bracketed numbers 
are the CI coefficients ol a given electronic configuration. Symbols 
H/L denote HOMO/LUMO orbitals. HF denotes the Hartree-Fock 
configuration. 

Peak E (eV) /12 Polarization Wave Function 











\HF) (0.8787) 










ft 


Y L-H ^ i)(0.1147) 


I 


4 1 


O.^J ( 00 


111 Ulcllic 


\H - 


> L + 2)(0.6566) 










\H - 


> L + 10) (0.3210) 










1 H" 


1 


L + 6) (0.2773) 










1 W 

— 


1 


i + 17) (0.1850) 




fills 


2.9911 


111 JJidlit; 


1 w 

— 


1 


i + 6) (0.5786) 










1 w 

— 


1 


L + 17) (0.3285) 










\H - 


3 ^ 


i)(0.2656) 










\H - 


2 ^ 


i + 11) (0.2544) 










\H - 


2 ^ 


i + 18) (0.2492) 


III 


6.639 


0.3735 


in-plane 


\H - 


1 


L + 3) (0.4496) 










\H - 


1 ^ 


i + 13) (0.4485) 










\H - 


4 ^ 


i) (0.4052) 










\H - 


> L4 


- 11) (0.2766) 










\H - 


> L4 


- 1)(0.2329) 










\H - 


> L + 18) (0.1905) 


IV 


7.311 


1.9267 


in-plane 


\H - 


^ L + 2)(0.3055) 










\H- 


2 ^ 


L + 11) (0.2892) 










\H - 


3 ^ 


L) (0.2834) 










\H - 


2 ^ 


i + 18) (0.2487) 










\H- 


2 ^ 


i + 1) (0.2029) 










\H - 


2 ^ 


i + 5)(0.1509) 


V 


7.842 


3.0105 


in-plane 


\H - 


> L + 2) (0.4233) 










\H - 




- 10) (0.3270) 










\H - 


2 


L -Ml) (0.2049) 










\H - 


^ L-h 20) (0.1946) 










\H - 


2 


L + l)(0.1594) 










\H - 


2 ^ 


i + 18)(0.1582) 



Note: "^GS does not correspond to any peak, instead it corresponds 
to the ground state wavelunction ol the isomer. 
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Table 7. Excitation energies, E, and many-particle wave functions of 
excited states corresponding to the peaks in the linear absorption spec- 
trum of square B4 (c/. Fig. 9), along with oscillator strength of transition 
(/12 = § ^(-£"2 — -El) X^i K'T^MilG) l'^), where, |m) denotes the excited state in 
question, \G}, the ground state, and dj is the i-th Cartesian component of the 
electric dipole operator. The polarization m-plane corresponds to the absorp- 
tion due to light polarized in the plane of isomer. In the wave function, the 
bracketed numbers are the CI coefficients of a given electronic configuration. 
Symbols H/L denote HOMO/LUMO orbitals. HF denotes the Hartree-Fock 
configuration. 

Peak E (eV) /12 Polarization Wave Function 



GS° 



I 4.879 4.3269 in-plane 



II 5.462 0.1792 in-plane 



III 6.418 0.2165 in-plane 



HF) (0.8682) 

H ^ L-H ^ L) (0.1765) 

H -2^ L-H -2^ L)(0.092Q) 

^L; J/ ^i + 1) (0.2946) 
H ^ L) (0.1909) 

H-2^L-H-3^ L)(0.1651) 

H ^ L;H ^ L + 1){0.5S93) 
H ^ L-H ^ L + 6) {0.2477) 

H ^ L) (0.5379) 

H ^ L-H ^ L + 1) (0.2922) 
H -2^ L;H -3^ L) (0.1744) 



IV 7.890 10.4628 in-plane |iy ^ L -M; - 1 L)(0.2351) 



Note: "GS does not correspond to any peak, instead it corresponds to the 
ground state wavefunction of the isomer. 
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Table 8. Excitation energies, E, and many-particle wave functions of 
excited states corresponding to the peaks in the linear absorption spec- 
trum of Hnear B4 (c/. Fig. 10), along with oscillator strength of transition 
(/12 = §^(i?2 — K'^Mi|G')P)j where, \m) denotes the excited state 

in question, \G), the ground state, and d.j is the i-th Cartesian component of 
the electric dipole operator. Longitudinal and transverse polarization corre- 
sponds to the absorption due to light polarized along and perpendicular to the 
molecular axis respectively. In the wave function, the bracketed numbers are 
the CI coefficients of a given electronic configuration. Symbols H/L denote 
HOMO/LUMO orbitals. HF denotes the Hartree-Fock configuration. 



Peak 


F, (pV) 
-C/ (_ev ; 


J 12 


J. VJlCXil IZjCXiLlKJlL 


Wave Function 






GS" 








\HF) (0.5636) 














\H L)(0.4737) 














L-H ^ 




- 1)(0.2291) 










\H - 1 ^ L; H ^ 


i4 


- 1)(0.2289) 


I 


5.363 


0.7433 


longitudinal 


\H ^ L;H -1^ 


i4 


-9) (0.3005) 












i4 


-9) (0.3005) 










\H ^ L; H - 1 ^ 


i4 


-8)(0.2702) 










\H ^ L; H - 1 ^ 


i4 


-8) (0.2702) 










\H ^ L; H - 1 ^ 


i4 


-4) (0.2439) 










\H ^ L; H - 1 ^ 


i4 


-4) (0.2439) 


11 


5.947 


31.8150 


longitudinal 


\H- 1 ^ L + 1)(0.2957) 










\H- 1 ^ L + 1)(0.2957) 


III 


7.352 


14.0546 


longitudinal 


\H ^ L;H -2^ 


L)(0.6477) 










\H-3^ L)(0.2966) 












\H- 1 ^ i + 1)(0.2091) 










\H- 1 ^ L + 1)(0.2091) 




7.380 


7.7738 


transverse 


\H- 1 ^ i + 13) (0.5253) 



Note: °GS does not correspond to any peak, instead it corresponds to the 
ground state wavefunction of the isomer. 
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Table 9. Excitation energies, E, and many-particle wave functions of excited 
states corresponding to the peaks in the linear absorption spectrum of dis- 
torted tetrahedron B4 (c/. Fig. 11), along with oscillator strength of transition 
(/i2 = §^(-£'2 — Ei)Y2i\{^\di\G)\'^), where, |m) denotes the excited state 
in question, |G), the ground state, and dj is the i-th Cartesian component of 
the electric dipole operator. The polarization x,y and z corresponds to the ab- 
sorption due to light polarized along x—,y— and 2— axis respectively. In the 
wave function, the bracketed numbers are the CI coefficients of a given electronic 
configuration. Symbols H/L denote HOMO/LUMO orbitals. HF denotes the 
Hartree-Fock configuration. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


GS" 








\HF) (0.6493) 










\H - 


s-L;Ji"-^L + 2) (0.4695) 










1 TJ 
\il — 


^ T . tj Q . r _L o\/'n 1 p;/17^ 


I 


1.000 


0.0517 


X 


\H - 


> L)(0.8127) 










\H - 


y L-H-1 ^ L + 2) (0.2364) 




1.111 


0.0714 


y 


1 n 

— 


T 1^ 12 — r L/ -f- If yU.Oy^'^ ) 










— 


i T ■ ft ^ T 4-l\('nc;i91"l 

7 1^^ 11 — 7 Lj -\- iJIU.OlZll 


II 


3.609 


0.2771 


X 


1 n 

— 


i T ■ ff i^r4-9\('n f^nQzL"! 

7 Iv, _o — 1 — ? 1^ -j- Afyu.Ouu^) 










1 w 

— 


7 Iv, _o — 1 — ? J^fyyJ.^OoO) 










1 n 






3.754 


0.3957 


y 


\H - 


>L + 1;//^L + 2)(0.2391) 


III 


4.48 


0.4818 


z 


1 T4 

— 




IV 


4.96 


0.5841 


z 


1 n 

— 


1 r\(r\ ofi7^^ 










1 n 




V 


5.92 


0.6954 


X 


\H - 


3 ^ L)(0.2191) 










\H - 


>L + 2;_H'-4^L)(0.2191) 


VI 


6.15 


0.6819 


y 


\H - 


> L;i?- 3 ^L + 1) (0.2526) 


VII 


6.508 


0.1217 


z 


\H - 


1 L-h 2) (0.7131) 










\H - 


1 L) (0.2104) 










\H - 


>L + 2;i?^L + 2) (0.2030) 


VIII 


6.858 


1.4336 


X 


\H - 


>L;iy-^L + 10) (0.6912) 










\H - 


^L + 2;ff^L + 10) (0.3415) 










\H - 


> L + 18;H ^ L) (0.2759) 




6.99 


1.1338 


y 


\H - 


> L;H ^ L + 5){0. 4169) 










\H - 


> L + 2;H ^ L + 5){0.3770) 










\H - 


> L + 3;H-1-^ L) (0.2403) 










\H - 


1 ^ L + 1) (0.2156) 


IX 


7.80 


0.5231 


z 


\H - 


1^L;H-1^L + 2)(0.3386) 










\H - 


^ L;i?- 4 ^L + 2) (0.3164) 










\H - 


>L;ii--4^L + 2) (0.2987) 










\H - 


1 ^ L;_ff - 1 ^ L)(0.2722) 










\H - 


^L + l;_ff-l^L + 3)(0.2174) 



Note: ° GS does not correspond to any peak, instead it corresponds to the ground 
state wavefunction of the isomer. 
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Table 10. Excitation energies, E, and many-particle wave functions of excited states 
corresponding to the peaks in the linear absorption spectrum of pentagon B5 (c/. Fig. 12), 
along with oscillator strength of transition (/12 = |^(i52 — ^-l) X^i K^-MjlC)!^)' where, 
|m) denotes the excited state in question, |G), the ground state, and di is the i-th Carte- 
sian component of the electric dipole operator. The polarization in-plane corresponds to 
the absorption due to light polarized in the plane of isomer. In the wave function, the 
bracketed numbers are the CI coefBcients of a given electronic configuration. Symbols 
H and L denote HOMO and LUMO orbitals respectively. HF denotes the Hartree-Fock 
configuration. 



Peak E (eV) 



/12 



Polarization Wave Function 











\H - 


1 - 


> L)(0.1906) 


I 


1.394 


0.1805 


_L to the plane 


\H - 


> L- 


f 1) (0.8536) 










\H - 


> L- 


f 11)(0.1779) 


II 


1.818 


0.0734 


in-plane 


\H - 


> i) (0.6524) 










\H - 


1 - 


> i?) (0.4900) 










\H - 


> L- 


f 3) (0.2218) 


III 


3.504 


0.6466 


in-plane 


\H - 


2 - 


-1- 1) (0.5804) 










\H - 


1 - 


> H)(0.5131) 










\H - 


> i) (0.3441) 


IV 


3.868 


0.2839 


in-plane 


\H - 


2 - 


^ L-h 1) (0.7968) 










\H - 


1 - 


> i?) (0.3416) 


V 


4.379 


0.1885 


in-plane 


\H ~ 


3 - 


> i?) (0.8185) 










\H - 


4 - 


> L) (0.1900) 


VI 


5.378 


1.2515 


in-plane 


\H - 


1 - 


> L) (0.7309) 










\H - 


4 - 


> L) (0.3583) 










\H~ 


1 - 


^ i + 3)(0.1994) 










\H - 


3 - 


> i?) (0.1808) 




5.576 


5.0600 


in-plane 


\H - 


2 - 


^ i-f 1) (0.4415) 










\H- 


3 - 


> L) (0.4045) 










\H - 


4 - 


> ff) (0.4000) 










\H - 


> L) (0.2399) 










\H- 


1 - 


^ i + 4) (0.2039) 



HF) (0.8541) 



Continued to Next Page 



Note: " GS does not correspond to any peak, instead it corresponds to the ground state 
wavefunction of the isomer. 
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Table 11. {Continued from last page) Excitation energies, E, 
and many-particle wave functions of excited states correspond- 
ing to the peaks in the linear absorption spectrum of pen- 
tagon B5 (c/. Fig. 12), along with oscillator strength of transition 
(/12 = §^(£2 — Ei)J2i where, |m) denotes the excited 

state in question, \G), the ground state, and di is the i-th Cartesian 
component of the electric dipole operator. The polarization in-plane cor- 
responds to the absorption due to light polarized in the plane of isomer. 
In the wave function, the bracketed numbers are the CI coefficients of 
a given electronic configuration. Symbols H and L denote HOMO and 
LUMO orbitals respectively. HF denotes the Hartree-Fock configura- 
tion. 



Peak 


E (eV) 


/12 


Polarization 


Wave Function 


Vii 


D.oUo 


4.yiDo 


in-planc 


\H - 


3 - 


> L) (0.4073) 










\H - 


1 - 


> L + 4) (0.3328) 










\H - 


1 - 


> L-h 9) (0.2952) 










\H — 


2 _ 












\H - 


4 - 


> L + 4) (0.2772) 










\H - 


4 - 


> L-H 9) (0.2531) 




6.528 


0.7269 


in-plane 


\H - 


1 - 


> L) (0.5596) 










\H - 


4 - 


> L) (0.3400) 










\H - 


2 - 


> L-h 2) (0.2761) 










\H - 


3 - 


> L + 4) (0.2667) 










\H - 


3 - 


> L + 9) (0.2598) 


VIII 


7.161 


1.9516 


in-plane 


\H - 


4 - 


> L + 4) (0.4305) 










\H - 


4 - 


> L + 9) (0.3859) 










\H ~ 


3 - 


> L) (0.3776) 










\H - 


1 - 


> L + 4) (0.2487) 




7.283 


0.5947 


in-plane 


\H - 


> L- 


f 6)(0.7102) 










\H - 


> L- 


f 9) (0.3359) 










\H - 


3 - 


> L + 4)(0.2091) 










\H ~ 


3 - 


> L-^ 9) (0.1854) 


IX 


7.702 


1.0247 


_L to the plane 


\H - 


> L- 


f 11)(0.7980) 










\H - 


> L- 


f 14) (0.2692) 




7.750 


0.8764 


in-plane 


\H - 


1 - 


> L + 3) (0.8332) 










\H - 


1 - 


> L-h 19) (0.1555) 
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Table 12. Excitation energies, E, and many-particle wave functions of excited 
states corresponding to the peaks in the linear absorption spectrum of distorted 
triangular-bipyramid B5 (c/. Fig. 13), along with oscillator strength of transition 
(/12 = §^(-^2 — where, \m) denotes the excited state in 

question, \G), the ground state, and dj is the i-th Cartesian component of the 
electric dipole operator. The polarization x,y and z corresponds to the absorption 
due to light polarized along x—,y— and z— axis respectively. In the wave function, 
the bracketed numbers are the CI coefficients of a given electronic configuration. 
Symbols H and L denote HOMO and LUMO orbitals respectively. HF denotes 
the Hartree-Fock configuration. 

Peak E (eV) /12 Polarization Wave Function 



VjrO 








\HF) (0.8615) 










\H- 


1 - 


^ //)(0.2165) 










1 Pt 

— 


1 

1 — 


^ ff- ff 1—4. r\(r\ T^yi"! 










— 




z — 


ff- ff 1 ^ r\(r\ iizL'^~i 


I 


1.774 


0.4172 


z 


\H - 


1 - 


-> L + 4) (0.8689) 


II 


3.161 


0.5817 


y 


\H - 


2 - 


-> //) (0.8734) 










\H- 


1 - 


^ L) (0.2303) 


III 


4.274 


0.3879 


z 


\H - 


1 - 


^ L + 2) (0.8666) 










\H- 


1 - 


-> L + 6) (0.1775) 


IV 


5.487 


1.4856 


y 


\H - 


1 - 


^ L) (0.7812) 










\H- 


2 - 


+ //)(0.2336) 










\H - 


2 - 


^ H-H -1^ L)(0.2332) 


V 


6.408 


0.8914 


z 


\H - 


4 - 


-> //) (0.7837) 










\H - 


2 - 


^ L + 2) (0.3357) 










\H - 


2 - 


^//;ff-l->L + l) (0.2274) 


VI 


7.519 


3.4471 


X 


\H - 


2 - 


^ H;H-1 ^ L + 2) (0.8133) 










\H - 


2 - 


^H-H-l ^ L + 6) (0.2669) 


VII 


7.744 


0.1349 


X 


\H ~ 


1 - 


^L;H-l -> L + 1) (0.6105) 










\H - 


1 - 


^ L + 4) (0.4305) 










\H- 


2 - 


-> L + 1) (0.2271) 



Note: GS does not correspond to any peak, instead it corresponds to the ground 
state wavefunction of the isomer. 
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